Analysis for SSFA project: Two factor model

Filtered strongly by < 20% NMP

Table of contents

  1. Used packages
  2. Global settings
  3. Load data
  4. Model specification
  5. Inference
    1. epLsar
    2. Asfc
    3. Smfc
    4. HAsfc9
  6. Summary

Used packages

Global settings

Output

Plotting

Computing

Load data

Show that everything is correct:

x1 indicates the software used, x2 indicates the treatment applied.

Model specification

Inference

epLsar

Verify model settings

Check prior choice

Prior choice is as intended: Broad over the data range.

Sampling

Save for later comparison

Check sampling

Posterior predictive distribution

Level plots

Posterior and contrasts

Verify model settings

Check prior choice

Sampling

Check sampling

Posterior predictive distribution

Compare prior and posterior for model parameters

Posterior and contrasts

Asfc

Verify model settings

Check prior choice

Prior choice is as intended: Broad over the data range.

Sampling

Check sampling

Posterior predictive distribution

Compare prior and posterior for model parameters

Posterior and contrasts

Smfc

Verify model settings

Check prior choice

Prior choice is as intended: Broad over the data range.

Sampling

Analysis stopped here because sampling did not converge. As the plot shows, some data points are very far away from the others, which would require the analysis to be based on more heavy-tailed distributions.

HAsfc9

Verify model settings

Check prior choice

Prior choice is as intended: Broad over the data range.

Sampling

Check sampling

Posterior predictive distribution

Compare prior and posterior for model parameters

Posterior and contrasts

Summary

Set the surface parameters for every treatment dataframe:

Show the treatment pairs and surface parameters where the softwares differ

Write out